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Abstract 

We study the generic NMSSM with the coupling XSH1H2 at the 
limit of perturbativity, A ss 0.7-^0.8, and a moderate stop mass, m t ~ < 
300 GeV. Respecting the LEP2 bounds and insisting on naturalness, 
CsJ , we delimit the allowed region of parameter space and we study the 

spectrum and phenomenology of the relatively light new particles. 
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1 Introduction 

Supersymmetry provides one of the most attractive solutions to the SM hier- 
archy problem. The missing evidence for a light Higgs, however, introduces 
the need for large loop corrections driven by the stop in order to keep the 
Minimal Supersymmetric Standard Model (MSSM) viable. The presence of 
a heavy stop requires some accidental cancellations to get the correct value 
of the Z boson mass, reintroducing a fine-tuning in parameters at the percent 
level. Moreover, this residual fine-tuning weakens the argument suggesting 
that supersymmetry should manifest itself around the weak scale, and thus 
be visible at the LHC. 

As it is well known, in the most straightforward extension of the MSSM, 
the Next to Minimal SM (NMSSM), there is an additional quartic term in the 
scalar potential that may help increasing the Higgs mass without the need 
to invoke large loop corrections (see [H El El HI [5] and references therein). 
The NMSSM contains indeed an extra gauge singlet chiral supermultiplet S, 
with the superpotential Yukawa interaction XSHiH 2 . In this framework, the 
mass of the SM-like Higgs, i.e. the one that couples to vector bosons before 
mixing, is given at tree level by: 

m\ = M 2 Z cos 2 2(3 + A V sin 2 2(3 . (1) 

It is evident from (CQ) that the value of the Higgs mass depends crucially on 
how big A is taken at the weak scale. In particular, if one does not want 
to spoil manifest perturbative unification, using the Renormalization Group 
equations one finds the upper bound A < 0.6 at low energies. Together with 
the complementary dependence on the anglqj (3 of the two terms in the r.h.s. 
of Eq. (p]), this makes again difficult to push rrih above the LEP2 bound 
without a sufficiently heavy stop. 

More recently it has been shown in [61 [TJ, IB] that if one assumes the 
presence of extra matter at intermediate energies, filling complete SU(5) 
supermultiplets, the RG evolution of A is slowed down. It is then possible 
to take A ~ 0.7 -j- 0.8 at the weak scale, obtaining a Higgs boson with a 
mass around 115 -r- 125 GeV, and a stop of moderate mass, consistently with 
perturbative unification. 

In [S] an explicit model was presented based on this framework, where 
a Peccei-Quinn symmetry in the superpotential was assumed, only weakly 
broken in the supersymmetry breaking terms. Other than solving the /i- 
problem, this allows to keep under control the number of new parameters. 
This is an interesting particular example which has, however, a restricted 
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range in parameter space where all experimental bounds are satisfied, without 
the need to tune some parameters at the 10 % level. 

In this work we consider the generic NMSSM with R-symmetry on the 
superpotential, while still keeping A ~ 0.7 -=- 0.8 at the weak scale and a 
moderate stop mass, m 4 ~ < 300 GeV. There is in this case the additional 
interaction kS 3 , that modifies the vacuum structure of the theory, opening 
up a wider region in parameter space. We consider the general breaking 
of the R-symmetry in the soft SUSY-breaking potential, including both the 
trilinear terms A^ and A\. Many recent works have considered the NMSSM 
with a spectrum of relatively light particles [9j QUI [TT] . We are not aware, 
however, of any work focussed on the generic NMSSM with A = 0.7 -j- 0.8, 
moderate tan (3 and a stop mass below 300 GeV. 

The paper is organized as follows. In Sec. [2] we present the model and 
we discuss the spectrum and the couplings of the Higgs scalars and of the 
Higgsinos. In Sec. [3]we discuss the naturalness of parameters and the bounds 
coming from negative searches at LEP2, and in particular their interplay in 
setting limits on the available parameter space. In Sec. H] we present the 
evaluation of the contribution to the S and T parameters from the extended 
Higgs sector. In Sec. [5] we discuss some of the experimental signatures of 
this model at the LHC. 

2 The Higgs-Higgsino sector 

In this Section we describe in detail the Higgs sector of the model, its pa- 
rameters, the mass spectrum and the main properties of the new particles. 

2.1 The scalar potential 

We consider the following superpotential: 

W = XSH X H 2 + ^S 3 . (2) 

The potential of Eq. (J2j) exhibits a continuous R-symmetry, that forbids 
any mass term in the superpotential, providing an elegant solution to the \x- 
problem. The S 3 term leads to relevant differences in phenomenology with re- 
spect to the PQ-symmetric limit considered in [8j. The soft supersymmetry- 
breaking potential, that includes R-symmetry breaking terms proportional 
to the couplings in W, is given by: 

V soft = mllH.l 2 + m 2 2 \H 2 \ 2 + m 2 s \S\ 2 + {XAxSH^ + \kA k S 3 + h.c.) . (3) 



Assuming CP conservation, the model contains seven free parameters: A, k, 
mf, m 2 ,, m|, A\, A^. When it exists, the CP-conserving symmetry-breaking 
vacuum respects the following minimization conditions: 

,9 9 m 2 — nil 9 2\v s ,, . . .. 

cos(2/3) sm(2/3j 

• <0R\ - 2Xv s (Ax + kv s ) 

Sm[P) ~ m? + m| + AV + 2A 2 ^' [) 

2k 2 vl + kA k v 2 s + v s (X 2 v 2 - Xkv 2 sin(2/3) + m|) - -AA A w 2 sin(2/3) = . (6) 

Using Eqs. (jU EJ E]) above, we can trade m\ and m 2 , for v and tan/3 = i^/ui- 
Moreover, three parameters are fixed by requiring the mass of the SM-like 
Higgs to be higher than the LEP bound. We choose A and tan/3 such that 
they allow rrih in Eq. (QQ) above 115 GeV. The value of A is determined by the 
request of perturbativity up to the GUT scale; once Xgut/^tt ~ 0.15 -j- 0.3 
is fixed, the RG running with extra matter (three 5 + 5 of SU(5)) gives 
the value at the weak scale A ~ 0.7 -r 0.8. Moreover, the request to have 
rrih > 115 GeV in Eq. (pQ) limits tan/3 in the interval 1.5 < tan/3 < 2.5, 
with a maximal value of rrih fc> r tan/3 ~ 2. The Higgs mass does not depend 
directly on k (before mixing). However, since the RG equations for A and k 
are coupled, it can be shown that a large k would reduce the maximal value 
of A. In order to avoid this effect, we have to choose k GUT /Air < 0.05; this 
correspond at the weak scale to k < 0.1. For a detailed discussion of these 
aspects see [S]. Throughout the rest of the paper, unless explicitly stated, 
we will consider A, k and tan/? fixed at the values that maximize the Higgs 
mass (A = 0.8, k = 0.1, tan/3 = 2). 

We are thus left with only three effective free parameters in the potential: 
m/g, A\, Ak- As it can be seen in Fig. [TJ for fixed A/, only a delimited region 
in the m| — A\ plane is allowed. This region shrinks for higher values of 
Ak- In Fig. [T] the bound on the right comes from global stability of the 
SU(2) x [/(l)-breaking minimum, while the ones on the left and top come 
from local stabilitya of the same minimum. The bound on the bottom is 
instead a consequence of the experimental limits on the chargino mass (see 
Sec. 



2.2 Spectrum 

The Higgs multiplets Hi,H2 and S appearing in the potential (J2J) contain 
seven bosonic physical degrees of freedom: two charged scalars, three neutral 
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Figure 1: Allowed regions (cyan shaded) in parameter space due to 
stability and chargino mass bounds, as a function of A\ and M = 

Sgn(m 2 s )J\m 2 s \, for A k = (left) and A k = 150 GeV (right). 



CP even and two neutral CP odd particles. Having denned 

H° = A={h + MTi), Hi = ±={h 2 + Z7T 2 ), S = -j={s + ITT,), (7) 

it is convenient to write down the scalar mass matrices in the basis where 
the eaten goldstone bosons decouple, given by: 



H = cos /3h 2 — sin/3/ii, h = cosj3hi + sin j3h 2 , s 



Pi = COS /?7T2 — Sin /?7Tl , P 2 



Tis 



(9) 



In this basis, among the CP-even fields only h has a trilinear coupling ghvv 
with the SM vector bosons. 

In the following we report the mass matrices of CP odd and even states, 
in the basis defined in Eqs. ((HI M>- We keep the explicit dependence of the 
matrix elements upon the singlet vev v s . Whenever it appears, v s has to be 
read as the numerical solution of Eq. (jHJ) as a function of the free parameters 
A\, A k and m 2 s . The squared mass matrix elements for the scalar degrees of 



freedom are, in the (H, h, s) basis of Eq. (jSJ): 

mix = ^r(A A + ^) + sin 2 (2/5)(m|-AV), 

m\ 2 = --sin(4/?)(m|-AV), 

m 2 3 = -\vcos(2/3)(A x + 2kv s ), 

m \ 2 = m \ cos 2 (2/5) + A V sin 2 (2/5) , 

m 2 3 = 2\ 2 v s v - Xv sm(2(3)(A x + 2kv s ), 

1 f 2 

m^ = -\A x sm(2(3)— + 4k 2 v 2 s + kA k v s . (10) 

The pseudoscalar squared mass matrix is, in the basis (Pi, P 2 ) of Eq. (jUJ): 



2Au s 



-(A A + fcy a ) Ai>(A A -2fa; a 



A^2 _ / sin(2£)' 

A ^ A«(A A -2Aw s ) ^l« 2 (f(A A + fct; s ) + 3AA;)-3M fc t; a 

'(H) 
The charged Higgs mass is given by: 

m 2 H± = m 2 w -\ 2 v 2 + X(A x + kv s )^-. (12) 

The resulting spectrum consists of: 

• two relatively light scalar particles, s\ and s 2 , the first lighter, the second 
heavier than the LEP2 bound on the SM Higgs of about 115 GeV; 

• a light pseudoscalar a\, in the 60 -f- 80 GeV, that would be massless in 
presence of a PQ symmetry; 

• three heavy, nearly degenerate states S3, a 2 , H^. 

The scalar particles spectrum is shown in Fig. [2] as a function of m|, for 
different values of Ak, and fixed A\ = 300 GeV. The effect of a moderate 
loop correction (mj = 300 GeV) is included. It has to be noticed that a 
positive A k increases the mass of the lightest scalar state Si, pushing down 
at the same time the light pseudoscalar a±. 

2.3 The Higgsino spectrum 

In the Higgsino sector, we assume that the gaugino mass parameters Mi and 
M 2 are large, so that gaugino-Higgsino mixing is negligible and the chargino 
mass can be kept above the LEP2 limit. In the basis 

Ni = ±=<fi x - H 2 ), N 2 = -^=(Hi+H 2 ), N 3 = S, (13) 
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Figure 2: Scalar spectrum as a function of M = S gn{m 2 s ) \ \rn? s \ for 
A k = 50 GeV (left) and A k = 150 GeV (right), and fixed A x = 300 
GeV. 



the Higgsino mass matrix is given by: 

/ n 



^§(cos/3-sin/?) \ 

fi -^(cos (3 + sin f3) 

\ ^(cos/?-sin/3) -^|(cos/5 + sin/3) ' " 2kv s , 







(14) 

The mass eigenstates Xi are found applying a rotation that diagonalizes the 
M. x matrix, such that: 



Xi = V ii N i VV 1 



1 VM X V T = M°. 



(15) 



The resulting mass spectrum consists of a light, mostly singlino, state and 
two heavier states, separated in mass by about 40 -r- 50 GeV (see Fig. [3]). 
The chargino mass is given (in the heavy gaugino approximation) by 

m x ± = \i. (16) 

The effective \x term appearing in (114j) and (TT6l) is given by 

fi = Xv s , (17) 

and thus no new parameters are introduced in the Higgsino sector. 

2.4 Couplings and decays 

In this section we discuss the properties of the Higgs sector particles rele- 
vant for phenomenology, such as production couplings and branching ratios. 
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Figure 3: Higgsino spectrum as a function of the chargino mass \x. 



A few characteristic properties have to be noticed, crucial to the collider 
phenomenology. The lightest scalar, si, acquires a trilinear coupling to vec- 
tors g Sl vv only through the mixing with the SM-like Higgs; the resulting 
suppression in production rates is shown in Fig. m where the squared nor- 
malized coupling £ 2 = gsusy/dsM * s plotted . Once produced, this particle 
will decay in light pseudoscalar pairs a\a\, whenever the channel is kinemat- 
ically allowed. Otherwise, it will decay predominantly in bb pairs, with the 
remaining part decaying in t + t~ with a relative weight given by 



BR{t+t- 
BR{bb) 



m b 



The light pseudoscalars a\ will in turn decay mostly in bb with the same BR 
ofEq. (USD. 

The scalar S2 has nearly full strength couplings to SM particles. Its 
production rates are therefore close to the SM ones. Moreover, the decay 
channel of S2 into neutralinos is always closed, due to the enhancing effect 
of the k term on the mass of the lightest neutralino. On the other hand, in 
wide regions of parameter space the decay into light pseudoscalar pairs aidi 
is dominant. Except for the region m| < —(80 GeV) 2 where the channel 
S2 — > S1S1 becomes also kinematically accessible, the fraction 1 — BR(s2 — ► 
aiOi) has SM-like decays. 

Another interesting property is that the two heavy neutralinos will decay 
with a sizable branching fraction in scalar particles (see Fig. [5]). This decays 
go through the singlino component of Xi ■ 
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Figure 4: Squared normalized couplings of s\. Coupling to vectors 
pairs £% y V (left) and top pairs ^ 2 # (right), plotted as a function of 

M = Sgn(m 2 s )J\m 2 s \ and A x , for fixed A k = 100 GeV. 



The branching fractions of the remaining heavy states are summarized in 
Figs. E] and El 

3 Parameter space: naturalness and LEP con- 
straints 

In this section we proceed further with the analysis of the parameter space, 
in order to check whether regions exist that are not yet excluded by the LEP2 
searches, consistently with a moderate fine-tuning on parameters. 



3.1 Naturalness 

Since one of the motivations to consider this model is the attempt to amelio- 
rate the naturalness problem of the MSSM, it is important to look carefully 
at all possible sources of fine-tuning and check if regions in parameter space 
exist where the residual fine-tuning is moderate q 

The first thing to check is the dependence of the weak scale v on the 
various parameters p«. This can be estimated by evaluating the logarithmic 



3 A similar analysis was performed in [16] : however no strict bounds where set on models 
with large A. 
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Figure 5: Typical branching ratios of heavy neutralinos xi (left) and 
X3 (right) as a function of M = Sgn{m 2 s )\J\m 2 s \. (A\ = 300 GeV, Ak = 
100 GeV). 



derivative of v 2 in Eq. 



<91ogt> 2 



d log p 2 



(19) 



taking into account the variation of (3 from Eq. (jSJ). The only relevant such 
dependence is found to be the one on A\. As it is shown in Fig. [TJ the request 
of a fine-tuning lower than ~ 10% disfavor values of tan/? smaller than 2 and 
A\ bigger than 400 -r 500 GeV. The other sizable source of fine-tuning comes 
from the one loop contribution to the running of m 2 s driven by A\. In the 
example of a PQ-symmetric potential studied in [8], it turned out that this 
was the main source of fine-tuning; the ratio 5m 2 s /m 2 s in fact was always 
bigger than 10, for a low mediation scale of A mess ~ 100 TeV. In the case of 
R-symmetric potential we are considering, there are instead sizable regions 
in the m 2 s — A\ plane, corresponding to negative m 2 s values, where this ratio 
takes smaller values. Fig. [8] shows the ratio 8m 2 s /m 2 s for A mess ~ 100 TeV. 



3.2 LEP limits 

The negative searches for supersymmetric particles performed by the LEP2 
experiments placed several lower bounds on the sparticle masses. Since the 
spectrum of our model contains some light particles, it has to be checked 
whether they are not already excluded by LEP, and what further restrictions 
on parameter space this implies. 

"* — > X1X2 and e + e~ — > 
Given the Higgsino 



In the Higgsino sector both the processes e + e 
— ► x + X~ have been searched for at LEP [12 
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Figure 6: Typical branching ratios of the heavy scalar S3 as a function 
of M = Sgn{m 2 s )J\m 2 \. (A x = 300 GeV,^ = 100 GeV). 



spectrum of our model (Fig. [3]), direct X1X2 production is not kinematically 
accessible. The limit on the chargino mass is m x ± > 103 GeV. Since m x ± = /i 
under the assumption of heavy gauginos, we obtain the bound in the bottom 
of m 2 s — A\ plane already shown in Fig. [TJ 

Another process that sets stringent restrictions on the parameters is 
e + e~ — > Zs\ [15] . As we discussed in sec. 12.41 once produced the light 
scalars s% decay either directly in bb pairs or in si — > a\a\ — > bbbb, depend- 
ing on available phase space. It is interesting to confront the regions in the 
m 2 s — A\ plane excluded by LEP with the ones compatible with a moderate 
fine-tuning on m 2 s . In fact the regions with big negative m 2 s where 5m 2 s /m 2 s is 
low, correspond to the s\ mass approaching the stability boundary m 2 3i = 0. 
One would thus naively expect that the same regions are already ruled out 
by LEP. It turns out, however, that for Ak in the range Ak ~ 50 -=-150 sizable 
region survive with moderate fine-tuning, as it is shown in Fig. [HI It has to 
be noticed that a low value of A^ is not unnatural since it is the breaking 
parameter associated to k, that as we already discussed has to be taken small 
by itself. 

Finally, the process e + e~ — > Z* — > 01S1 is also constrained by LEP, but, 
due to the extreme weakness of the gz Siai coupling, it does not set further 
bounds on parameter space. 
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Figure 7: A^ A as a function of A\ for Ak = 100 GeV and m| 
— (60 GeV ) 2 . This quantity is almost independent on m| and Ak- 



4 Electroweak Precision Tests 



Besides the direct experimental constraints discussed in Sec. 13.21 we should 
also consider the contributions that arise from the extended Higgs-Higgsino 
sector to the electroweak precision parameters, S and T. The presence of a 
sizable Yukawa coupling A ~ 0.8 may in fact induce significant effects. 
The contributions to S and T from scalar particles are given by: 



T scal = J2Zvv Si 3{ A (m Si ,m z ) - A(m Si ,m w )) + ^2^ H±s .A(m H ±,m Si ) 



i=l 

2 



i=l 



3 2 



+ J2twH± aj A ( m H±,m Si ) -^2^2^ 2 Zs . a .A(m aj ,m Si ), (20) 

i=\ i=i 



3=i 



^scal = ^&Vsfi( F ( m s^rn z )+m 2 z G{m Si1 m w )) 



3 2 



+ SX^^K^mJ - F(m H ±,m H ±) 
i=i j=i 



(21) 
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Figure 8: 5m 2 s /m 2 s as a function of M 
(A k = 100 GeV). 



Sgn(m 2 s )^\ 



mV\ and A\ 



iwsi = U i2 


iwH± — Un 


iz &l a x = U ix cos a 


iz Sl a 2 = Un sin a 


iwH± ai = cos a 


iwH±a 2 = sin a 



Table 1: Reduced scalar- vectors couplings (£ = gsi/SY /9sm) 



where s,- 



are the scalar (pseudoscalar) mass eigenstates, the reduced 



couplings £ are defined in Tab. HJ while the complete expressions for the 

loop functions A, F, G can be found in App. \M 

The coefficients Uij in Tab. Q] are the elements of the matrix that rotates 

from the basis (H, h, s) of Eq. (jSJ), where only h couples to vectors, to the 

mass eigenstates basis. The rotation angle a is defined in an analogous way 

as: 



ax = — sma7i s + cosa(cos/?7r 2 — sin /3ttx) 
a>2 = cosa7r s + sina(cos/?7T2 — sin/?7Ti). 



(22) 
(23) 
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Figure 9: Contours of the ratio 5m 2 s /m 2 s in the region not excluded by 
LEP, as a function of M = S gn{m 2 s ) \ \m 2 s \ and A\, for A^ = 50 GeV 
(left) and A k = 150 GeV (right). 



The Higgsino contributions to S and T are given by 

v3 



-m* 



Tffiggsinos = J2 1 ( V £ ) 2 ^' m *) + ( V 2 T i) 2 M^ 

- \ HljJ V S V S + VM^Mm, -mj), (24) 

\ E-^i^J + ^iT^) 2 FK, -mj) - F(/i,/i), (25) 



S 



Higgsinos 



where [i is the chargino mass, the Higgsino/neutralino rotation matrix V has 
been defined in f fl~5"]) . and the complete expressions for the loop functions 
A, F can be found in App. [A] 

In order to compare our results with experimental constraints, we want 
to place the contributions from the model in the S-T plane. This is done in 
two steps: 

1. The contributions of a SM Higgs 



7ffig gs (>ft) = -3 [A(m h ,m w ) - A(m h ,m z )} , 
Smggs{m h ) = F(m h ,m z ) +m 2 z G(m h ,m z ) 



(26) 
(27) 



are subtracted from the SM S-T values of p2] , for a reference value of 
the Higgs mass; 

2. The contributions ( 1201 I2T1 1241 1251) from our model are then added. 
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Fig. [TU] shows the result as a function of A\ and tan/3, for A = 0.8. As can 
be clearly seen the S and T values of our model are fully compatible with 
the experimental contours, for all relevant values of the free parameters (the 
dependence of S and T on Ak and m| is indeed rather mild). 
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Figure 10: S and T values as a function of A\ (in the range 200 GeV - 
600 GeV) and tan /3, compared to the SM result and experimental limit 
(A k = 100 GeV, m§ = -(60 GeV) 2 ). 



5 LHC Phenomenology 

In this section we discuss some of the experimental signatures that should 
be looked for at the LHC, if this model is realized in Nature. We will mainly 
focus on the most promising channels for SUSY discovery at low luminosity. 
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A preliminary remark is the following. As indicated in Sec. 13.14 natu- 
ralness considerations prefer (although only logarithmically) a relatively low 
mediation scale A mess , and in turn a low scale of supersymmetry breaking 
y/F. In this case a light gravitino is present, that takes the role of the LSP, 
while the lightest neutralino is the NLSP. A crucial property for phenomenol- 
ogy is thus the lifetime of the Xi- If this particle is long lived, it will decay 
to gravitino well outside the detector and can thus be considered in the ex- 
perimental signatures as the effective LSP; otherwise, if it decays inside the 
detector, the phenomenology will appear quite different. 

For an Higgsino-like neutralino, as we assumed throughout this paper, 
the only possible decays are Xi ~^ 4>G, where G is the gravitino and </> a 
generic scalar. In particular, given the scalar spectrum of our model, the only 
kinematically allowed 2-body decay is the one involving a light pseudoscalar 
a\. The decay width is given by 



4 X1B1 ™xi A m li 



32tt (VF) 4 V m l 



XI 



(28 



with the coupling g^ a given bya 

g 2 ^ a = | cosa(V A n sin/5 — V12 cos/5) + V13 sina| 2 . (29) 

The neutralino decay length can be expressed as 

1.97 x 10- f & XV. , ™^- (Jf\ ( _ ^ - 



g% \ml J VlOOGeV/ I 100 TeV / V ml 

(30) 

where E is the neutralino energy in the lab frame. Depending on the SUSY 
breaking scale y/F, Eq. ( |30l) leads to a decay inside or outside the detector. 
The critical value for y/F is a few hundreds TeV. 

In the following we will refer to the case where the Xi lifetime is long 
enough such that it decays outside the detector and can be considered like 
the effective LSP. Otherwise, the experimental signatures will include four 
additional fo-jets (originated in two displaced secondary vertexes) coming 
from the decay chain xi ~> a \G —* bbG. These additional jets might allow 
for a more efficient signal/background discrimination. However the discussion 
of such complicated topologies is beyond the scope of this work. 

The interesting signatures for SUSY discovery depend on the full spec- 
trum of supersymmetric particles, including gluino, squarks and sleptons. 



1 Mixings are denned in Eqs. (J2UHS]). 
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In order to ease possible problems in the flavor sector, and to simplify the 
discussion, we will assume that all those particles are heavy (in the 1-2 TeV 
range), except for the ones that cannot be taken heavy without introducing 
additional fine-tuning [IS] . As already pointed out, the strongest constraint 
is the one on the stop mass. If we want to avoid big corrections to the soft 
mass parameters, we have to assume m t ~ < 300 GeV. Furthermore, if we look 
at the coupled RG equations for the stop and gluino, to avoid accidental 
cancellations we have to take m- g < 2 -j- 2.5 m 4 ~. 

If this interesting limit is realized, the most promising channel for SUSY 
discovery is the gluino pair production. For a gluino mass in the range 
400 -T- 600 GeV the production cross section at the LHC varies between 200 
and 10 pb [H]. The produced gluinos then decay with 100% BR in top-stop 
pairs, since the stop is the only light squark . The stop decay pattern depends 
on the details of the stop mixing matrix. As it is shown in Fig. [Tl] , if the 
lightest stop is mostly Ir or maximally mixed, it will decay mostly in b\ + ■, 
with the chargino further decaying as x + ~^ W/W*Xi- Otherwise, if the t\ 
is mostly ti it will decay in t\i with ~ 100% branching ratio. The final 
state thus consists in four 6-jets, missing energy and some combination of 
jets and leptons originating by the intermediate state W bosons. The main 
background sources for this process are ti, W/Z+jets and 66+jets. In order 
to reduce the background cross-sections the particular structure of the final 
state can be exploited, by imposing cuts on the appropriate variables, such 
as the number of 6-tagged jets, B]t and M e ff. These selection cuts should 
allow to extract signal from background with a luminosity of the order of 1 
fb _1 (this preliminary conclusion is supported by the study of similar decay 
chains in [2H |22] , where however a different mass range was considered) . A 
careful analysis of this process would be of much interest, in order to assess 
on firm grounds the SUSY discovery potential for this model. 

Obviously, a crucial experimental issue is the search for the light scalars, 
Si and S2- As was already pointed out in Sec. 12.41 S2 has in a wide regions of 
parameters the dominant decay S2 — > OiOi — » bbbb. This channel is likely to 
require a high luminosity, even if the gluon fusion cross section is unaffected 
by the mixing, because of the overwhelming SM backgrounds. The discovery 
potential for this mode, however, could be enhanced if the associated W/Zh 
production is considered, as it is pointed out in [20]. Other channels worth 
to be considered, as usual for light scalars, are S2 — > ZZ* /WW* and the loop 
induced decay S2 — > 77. Even with a small branching ratio, these channels 
lead to very clean experimental signatures. The observation of the light scalar 
Si seems to be a challenging task, because due to the suppression of the g si vv 
and g Slt t couplings (see Fig. 0j) there is a sizable reduction of the production 
cross sections as well as the decay rates in the channels mentioned above. A 
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Figure 11: Branching ratio of t\ — » 6% + f° r *i ~ ^R (left) and ti ~ 
(l/y/2)(ii, + in) (right), as a function of M = Sgn{m 2 s 
(A k = 100 GeV). 
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detailed study is needed to assess the discovery potential for both the s x and 
s 2 particles. 

Finally, another interesting process is the direct neutralino-chargino pro- 
duction pp — > X2X ± ■ The expected production rate for this process is reported 
in Fig. [T2l Without entering into details, it can be roughly estimated that 
this process would require a luminosity of at least 10 fb _1 in order to be 
detected in the trilepton channel. 



6 Conclusions 

In this work we studied the generic NMSSM with an R-symmetry invariant 
superpotential, in the regime where the XSHiH 2 coupling is taken at the 
limit of perturbative unification. Under the assumption of extra matter at an 
intermediate scale, the A value at the weak scale can be as big as A ~ 0.7-^0.8, 
allowing to raise the Higgs mass without the need of big radiative corrections. 
We insist in fact on naturalness, and we thus consider a stop of moderate 
mass, mi < 300 GeV. 

The R-symmetry in the superpotential forbids explicit mass terms, solv- 
ing the /i-problem, but allows for the interaction kS 3 besides the A term. 
We consider the soft supersymmetry-breaking potential that breaks the R- 
symmetry in a general way, including trilinear terms associated to both the 
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A and k supersymmetric couplings. Requiring to have the Higgs boson mass 
in its maximal range, 115 -r- 125 GeV, we are left with only three effective 
free parameters: A\, Ak and m|. 

The Higgs boson spectrum contains two light scalars, si and s 2 , one 
above and one below the critical value of 115 GeV, a light pseudoscalar ai, 
and three heavy states (the scalar s 3 , the pseudoscalar a 2 and two charged 
particles H^). The two light scalars, s± and s 2 , share the SM coupling to 
vector bosons pairs gtvv- Both si and s 2 prefer to decay in pseudoscalar 
pairs aiai, whenever this channel is kinematically open. In the neutralino- 
chargino sector we assume heavy gaugino masses. In this limit the light 
states, three neutralinos and one chargino, have only Higgsino components. 
For the remaining supersymmetric particles, we consider the interesting limit 
where all of the spectrum is heavy except for the states that have to be kept 
light for naturalness reasons. We thus assume, as already pointed out, the 
lightest stop to have a moderate mass (m^ < 300 GeV) and a gluino not 
exceeding twice of the ti mass, while all the remaining squarks and leptons 
are in the 1 -¥■ 2 TeV range. 
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We delimit the allowed region of parameter space where all the LEP2 
bounds are respected, insisting on the request to keep the fine-tuning mod- 
erate. The most stringent experimental bound is the one on the light scalar 
sx- However, since sx has a suppressed coupling to vector pairs, there are big 
regions in parameter space where it would not have been detected at LEP2 
even if it is lighter than 115 GeV. On the other hand the naturalness anal- 
ysis shows that wide regions in parameter space exist where the fine-tuning 
can be kept below the ten percent level, for a relatively low supersymmetry 
breaking scale of some hundreds of TeV. For Ak in the range 50 < Ak < 150 
GeV, there are sizable regions in the m 2 s — A\ plane where the experimental 
bounds are respected consistently with a moderate fine-tuning. 

We also consider the contribution to the Electroweak Precision param- 
eters S and T from the extended Higgs-Higgsino sector of the model, since 
A ~ 0.8 could introduce a sizable effect. We find that these contributions are 
perfectly consistent with the experimental data. 

The LHC phenomenology of the model depends crucially on the lightest 
neutralino properties. Naturalness suggests a low supersymmetry breaking 
scale y/F, so that the gravitino is the LSP, while the Xi is the next-to-LSP. 
The width of the xi ~ > a iG decay determines the lifetime and the decay 
length of the NLSP. For \JF greater than a few hundreds TeV, the xi wm 
decay outside the detector, and can be thus considered like an effective LSP 
in the experimental analysis. In this case we expect SUSY discovery in the 
gluino pair production process with ~ 1 fb _1 of luminosity or less, depending 
on the gluino mass. Other interesting processes at higher luminosity are the 
production of the light scalars and Higgsinos. A detailed analysis of these 
signatures would be of great interest to assess the LHC discovery potential 
for this model. 
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A One loop contributions to S and T 

We report here for the reader's convenience the expressions [2HI |2"1] of the one 
loop function that appear in the contributions to S and T from new particles. 
For a boson loop with internal masses mi,m2, the functions A, F, G appearing 
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in Eqs. ( |2"U} |2"T]) are given by: 
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For a fermion loop with internal masses mi,m2, the functions A, F appearing 
in Eqs. ( 1241 1251) are given by: 
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